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Abstract 

We discuss some formal and fundamental aspects related with the replacement of 
integral dispersion relations by derivative forms, and their practical uses in high 
energy elastic hadron scattering, in particular pp and pp scattering. Starting with 
integral relations with one subtraction and considering parametrizations for the to- 
tal cross sections belonging to the class of entire functions in the logarithm of the 
energy, a series of results is deduced and our main conclusions are the following: (1) 
except for the subtraction constant, the derivative forms do not depend on any addi- 
tional free parameter; (2) the only approximation in going from integral to derivative 
relations (at high energies) concerns to assume as zero the lower limit in the inte- 
gral form; (3) the previous approximation and the subtraction constant affect the 
fit results at both low and high energies and therefore, the subtraction constant can 
not be disregarded; (4) from a practical point of view, for single-pole Pomeron and 
secondary reggeons parametrizations and center-of-mass energies above 5 GeV, the 
derivative relations with the subtraction constant as a free fit parameter are com- 
pletely equivalent to the integral forms with finite (non-zero) lower limit. A detailed 
review on the conditions of validity and assumptions related with the replacement 
of integral by derivative relations is also presented and discussed. 

Key words: elastic hadron scattering, dispersion relations, high energies 
PACS: 13.85.Dz, 13.85.Lg, 13.85.-t 



To appear in Nuclear Physics A 



Preprint submitted to Elsevier Science 



1 February 2008 



Contents 

1. Introduction 

2. Historical summary and main points 

3. Derivative dispersion relations 

3.1 Formal derivation 

3.2 Discussion 

4. Tests with Pomeron-Reggeon parametrizations 

4.1 Basic comments 

4.2 Degenerate meson trajectories 

4.2.1 Analytical results 

4.2.2 Global fits 

4.3 Non-degenerate meson trajectories 

4.4 Discussion 

5. Conclusions and final remarks 



1 Introduction 

Elastic hadron-hadron scattering, the simplest soft diffractive process, con- 
stitutes one of the hardest challenges in high-energy physics. Despite all the 
success of QCD concerning hard and also semi-hard processes, the large dis- 
tances involved in the elastic sector (and soft processes in general) demand 
a nonperturbative approach and, presently, we do not know how to calculate 
an elastic scattering amplitude in a purely nonperturbative QCD context. 
At this stage, phenomenology is very important, but must be based on the 
constraints imposed on the scattering amplitudes by some rigorous theorems, 
deduced from general principles of the underlying local quantum field theory, 
namely, Lorentz Invariance, Unitarity, Analyticity and Crossing. 

Dispersion relations play an important role in several areas of Physics, both 
as a practical tool and a formal theoretical result. In special, for particle- 
particle and particle-antiparticle interactions, they are consequences of the 
principles of Analyticity and Crossing [1]. In this context, they correlate real 
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and imaginary parts of crossing even (+) and odd (— ) amplitudes, which in 
turn are expressed in terms of the scattering amphtudes for a given process 
and its crossed channel, for example, a + b and a + b: 

F,,^F+ + F_ F,i^F+-F_. (1) 



Among the physical quantities that characterize high-energy elastic hadron- 
hadron scattering, the total cross section (optical theorem) and the p param- 
eter (related with the phase of the amplitude) are just expressed in terms of 
forward real and imaginary parts of the amplitude [2], 

, , Im F(s,t = 0) 
P^'^ = ImF(.,t^O) ' 



where s and k are the center-of-mass energy squared and the momentum, re- 
spectively, and t is the four-momentum transfer squared. Therefore, a natural 
and well founded framework for investigating the behaviors of cxtot and p, as 
function of the energy, is by means of dispersion relations. On one hand, from 
analytical parametrizations for (Jtot and fits to the corresponding experimental 
data, the p parameter may be determined (analytically or numerically) and 
compared to the experimental data and/or extrapolated to higher energies. 
On the other hand, cxtot and p may be analytically connected through the dis- 
persion relations and, therefore, both may be determined from global fits to 
the experimental data (that is, simultaneous fits to (Jtot and p data), improving 
the statistical ensemble in terms of degree of freedom (F) . 

For particles and antiparticles, the pp and pp scattering correspond to the 
highest energy interval with available data, reaching y/s ~ 63 GeV for pp 
(CERN ISR) and ^/s ~ 2 TeV for pp (Fermilab Tevatron). From experiments 
that are being conducted at the Brokhaven RHIC it is expected data on pp 
scattering at s/s: 50 - 500 GeV, and in the near future the CERN LHC will 
provide data on pp scattering at 16 TeV. It is highly expected that, in short 
term, these novel experiments will allow crucial tests for several phenomeno- 
logical models and, among them, the analytical models, based on the use of 
dispersion relation techniques, certainly play a central role. 

The use of dispersion relations in the investigation of scattering amplitudes 
may be traced back to the end of fifties, when they were introduced in the form 
of integral relations. Despite the important results that have been obtained 
since then, one limitation of the integral forms is their non-local character: in 
order to obtain the real part of the amplitude, the imaginary part must be 
known for all values of the energy. Moreover, the class of functions that allows 



3 



analytical integration is limited. By the end of the sixties and beginning of 
the seventies, there appeared the dispersion relations in a differential form, 
and they provided new insights in the dispersion relation techniques. However, 
as we shall discuss, there still remain some questions related with both the 
formal replacement of integral by the derivative forms and their practical uses 
in high energies, as, for example, the effects of the approximations considered, 
the correct expression of derivative form (depending on the class of functions 
involved) and the role of the subtraction constant. We understand that these 
aspects must be made clear before any practical use of the derivative relations. 

In this work we critically discuss the situation concerning the derivative re- 
lations and present some answers for the above mentioned questions. The 
manuscript is organized as follows. In Section 2 we briefly review some his- 
torical facts and results about the integral and derivative relations, stressing 
the problems we are interested in. In Section 3 we first demonstrate that, for 
the class of functions that are entire in the logarithm of the energy, except 
for the subtraction constant, the derivative relations do not depend on any 
additional free parameter. We then review the main results in the literature, 
obtained through different approaches, in connection with the corresponding 
assumptions and classes of functions involved. In Section 4, making use of 
a pomeron-reggeon parametrization, we analyze the practical effects of the 
approximations involved in the replacement of the integral by the derivative 
forms. In particular, we show that, for this class of functions, the derivative 
form with the subtraction constant as a free fit parameter is equivalent to the 
integral form for ^/s above 5 GeV. The conclusions and some critical remarks 
are the contents of Section 5. 



2 Historical summciry and main points 

In this Section we briefly review some results related with the replacement of 
integral dispersion relations (IDR) by derivative dispersion relations (DDR), 
stressing the points that we shall discuss in the next sections. We are interested 
in the high-energy region, specifically ^/s > 5 GeV and, as commented in 
our introduction, the practical tests will be performed with the pp and pp 
experimental data on atot and p. For these scatterings, the experimental data 
above \/i ~ 20 GeV indicate a slow increase of cxtot, roughly as ~ In^ s and that 
afot — crfot ~ as s — > oo. These conditions allow the use of integral dispersion 
relations with only one subtraction, and in the standard form they are usually 
expressed in terms of the energy E and momentum p of the incoming proton 
in the laboratory system [2,3]. Since we shall be interested in the region ^/s > 
5 GeV, we can approximate s = 2m{E -|- m) ~ 2mE with an error < 5%, 
which decreases as the energy increases. In this case the standard IDR, with 
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poles removed, are given by 



ReF+(s) 



1 



ImF+(s'), 



(4) 



TT J 

so 





ReF_(s) 




ImF_(s'), 



(5) 



where K is the subtraction constant, and for andpp scattering, sq = 2rn? ~ 
1.8 GeV^. It should be noted that for the even amplitude the relations with 
one and two subtractions are equal and for the odd case the relations without 
subtraction and with one subtraction are equal. 

Through Eqs. (1) to (5) the physical quantities Utot and p may be simultane- 
ously investigated. By means of IDR, important results have been obtained 
since the beginning of the seventies, as for example in the works by Bourrely 
and Fischer [4], Amaldi et al. [5], Block and Cahn [2], Kluit and Timmer- 
mans [6], Augicr et al. (UA 4/2 Collaboration) [7], Kang, Valin and White 
[8], Bertini et al. [9] and many others. However, as commented in our intro- 
duction, the inconveniences with the IDR concern their non-local character 
and the limited number of functions that allows analytical integration, and 
therefore, error propagation from the fit parameters. On the other hand, un- 
der some conditions, the above integral forms may be replaced by quasi-local 
ones, expressed in a derivative form and called derivative dispersion (or ana- 
lyticity) relations. The first result in that direction appeared indicated in the 
works by Gribov and Migdal [10], and afterwards, the DDR were treated in 
more detail by Bronzan [11], Jackson [12], Bronzan, Kane and Sukhatme [13]. 
In particular, the forms deduced by the last authors are given by 



where a is a real parameter. Soon after, based on the Sommerfeld-Watson- 
Regge representation and other assumptions, Kang and Nicolescu introduced 
the following expression for the derivative relations [14] : 
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(9) 



which, according to the authors, reduces to the Bronzan, Kane and Sukhatme 
result for a decreasing power form of Im[F_ (s)/s] [14]. In the last years, these 
equations have been extensively used in the detailed analysis on analytical 
models by Cudell et al. [15] and by the COMPETE Collaboration [16]. 

The DDR, in the form introduced by Bronzan, Kane and Sukhatme (BKS), 
have been criticized by several authors [17,18,19,20] and analyzed and dis- 
cussed in detail by Fischer, Kolaf [21] and Vrkoc [22]. Despite all these works, 
wc understand that two aspects related with the replacement of the IDR, Eqs. 
(4) and (5), by the DDR, Eqs. (6), (7), or (8), (9) and with their practical 
uses are yet unclear in the literature. 

One aspect concerns the origin, role and range of the exponent a in the BKS 
relations, which does not appear in the formulas by Kang and Nicolescu. For 
example, in the paper by BKS this parameter is allowed to have any real 
value [13], but in References [18,23] (see also [21]) it must lie in the interval 
< a < 2. According to some authors, "the choice of a different from 1 has no 
practical advantage" [21] and other authors have used it as a free fit param- 
eter [13,24], which demands some physical interpretation. The calculation by 
Eichmann and Dronkers, in 1974, lead to DDR that do not depend on a [17], 
but this parameter appears in a recent review by Kolaf and Fischer [25]. To 
our knowledge, with the exception of the above quoted considerations of a as 
a free fit parameter, all the other practical uses of the BKS relations are made 
by assuming a = 1. However, a — 1 leads to a singularity in the BKS expan- 
sion of the odd amplitude, since it contains the term sec^(7ra/2) [13]. As we 
shall discuss in detail in Section 3.2, some of these results and statements are 
not contradictories, because they refer either to different energy regimes (finite 
or asymptotic s) or to different classes of functions, for example, entire func- 
tions in the logarithm of the energy, or functions satisfying general principles 
from axiomatic field theory, etc... Moreover, some results have been obtained 
in different historical contexts in terms of the highest energies reached in ex- 
periments, leading to different mathematical assumptions on the asymptotic 
conditions. 
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Another aspect concerns the subtraction constant in the singly subtracted 
IDR, Eq. (4), which is the starting formula in the work by Bronzan. Kane and 
Sukhatme. That constant does not appear in all the above derivative forms 
and, to our knowledge, neither in almost all their practical uses, including the 
detailed works by the COMPETE collaboration. Yet, to our knowledge, the 
only exception concerns the recent works [26,27,28], where we have performed 
analysis of the experimental data available, including cosmic-ray estimations 
for the total cross sections, by using DDR with the subtraction constant as 
a free fit parameter and have shown that it affects the fit results at both low 
and high energies; therefore, the subtraction constant can not be disregarded. 
We have also shown that, for functions that are entire in the logarithm of the 
energy (Taylor expansion), the derivative forms do not depend on the param- 
eter a [28]. More recently, Cudell, Martynov and Selyugin have introduced a 
new representation for the DDR, intended for low energies as well, and they 
also refer to the possible role of the subtraction constant in the derivative form 



Based on the above considerations, in the next sections we investigate some 
formal and practical aspects related with the replacement of the IDR by the 
DDR. Specifically we shall demonstrate that, for entire functions in the loga- 
rithm of the energy, the derivative forms do not depend on the a parameter 
and shall discuss both the approximations involved and the important prac- 
tical and formal role of the subtraction constant. Some of these results have 
already been presented elsewhere [28]. 



3 Derivative dispersion relations 

In this section we first show that the subtraction constant is preserved when 
the IDR are replaced by the DDR and that, for functions entire in the log- 
arithm of the energy, the derivative forms do not depend on any additional 
free parameter. Next, based on the formulas displayed, we present a detailed 
discussion on the conditions of validity of the different results mentioned in 
the previous section, in connection with the assumptions involved. 

3.1 Formal Derivation 

Let us consider the even amplitude, Eq. (4). By defining s' = e^', s = e^ and 
g{^') = ImF+(e^ )/e^ , we express 
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Assuming that g is an analytical function of its argument^ we perform the 
expansion 



9{0 = t ^9(0 

n=0 '-^'S 



nl 



and then, the integration term by term in the above formula. At the high 
energy hmit, we consider the essential approximation sq — 2m^ — > 0, so that 
In So — oo. With these conditions, we obtain 



sinh(e - 



Now, defining y = C — the above formula may be put in the form 
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and, therefore, 
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where the series expansion is imphcit in the tangent operator. With analogous 
procedure for the odd relation we obtain 



ReF_(s) = tan 
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We see that, with the exception of the subtraction constant, these formulas 
do not depend on any additional free parameter, or, they correspond to a 
"particular" case of the BKS relations for a = 1. We stress that the same 
result is obtained if, following BKS, the derivation begins with an integration 
by parts without free parameter, as can be easily verified. We shall return to 
this point in Section 3.2. 

In practical uses the trigonometric operators are expressed by the correspond- 
ing series. With the substitution in the odd case 
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and, by expanding the series around the origin, we obtain 
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which is formally equivalent to the results introduced by Kang and Nicolescu, 
Eqs. (8) and (9). except for the presence of the subtraction constant. These 
expansions have been used in the Ref . [26] , that includes power and logarithmic 
dependences in the parametrizations. 

We have just demonstrated three formal and important results: (1) the sub- 
traction constant is preserved when the IDR are replaced by DDR and, there- 
fore, in principle, can not be disregarded; (2) except for the subtraction con- 
stant, the DDR with entire functions in the logarithm of the energy, do not 
depend on any additional free parameter; (3) the only approximation involved 
in the replacement concerns the lower limit in the IDR, namely, Sq = 2m^ 0, 
which represents a high-energy approximation. We shall return to these points 
in Section 4, when treating practical uses of the DDR. 

3.2 Discussion 

With focus on the DDR as expressed by Eqs. (6) to (13), let us discuss here the 
different results mentioned in Section 2, in connection with the corresponding 
conditions of validity and assumptions involved. We shall roughly follow a 
chronological approach (see also [25]). 

The first derivative form appeared in the beginning of 1968 in the context 
of the Regge Theory. By investigating the cuts associated with the Pomeron 
in the high energy limit [s — > oo) and at low momentum transfers, Gribov 
and Migdal introduced a derivative representation for the Watson-Sommerfeld 
integral, which, mathematically, corresponds to the first term in the expansion 
(12) (without the subtraction constant), namely 



In the beginning of the seventies, Bronzan [11], followed by Jackson [12], in- 
troduced the first form involving the tangent operator, Eqs. (10) and (11), 
independent of the parameter a. The arguments by Jackson were based on 
the formal representation of the Taylor series by an exponential operator. 



an on the high energy limit, in order to avoid branch points and cuts (see [12] 
for details). Therefore, the basic assumptions concern entire functions and the 
high energies. 
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The first direct connection between IDR and DDR appeared in 1974 in the 
work by BKS [13]. Starting from the singly subtracted IDR, Eq. (4) (neglecting 
the subtraction constant), the authors mention an integration by parts, which 
leads to an expression containing the parameter a. That may be obtained 
by multiplying and dividing Eq. (4) by and then integrating by parts. 
In addition to the high-energy approximation, represented by the lower limit 
So 0. a Taylor expansion of Im / is also assumed, so that the series may 
be integrated term by term (uniform convergence), leading to Eqs. (6) and (7). 
We should note that the expansion considered by BKS in [13] is around the 
point T:{a — l)/2 for the even amplitude and T:a/2 in the odd case, which 
leads to the previously mentioned singularity for a = 1. Differently, Eqs. (8), 
(9) and (12), (13) and all the other cases treated here and discussed in what 
follows, refer to expansions around the origin. 

In the same year, Kang and Nicolescu introduced the series expansion (8) 
and (9), which, according to the authors, can be derived "from the analytic- 
ity relation given by the Sommerf eld- Watson- Regge representation" and are 
asymptotic relations [14]. The point here was to treat the possibility that 
the amplitude could include not only simple poles in the complex angular 
momentum plane but, in particular, logarithm dependences in the odd am- 
plitude, representing the Odderon [30]. The association of this concept with a 
complicated singularity at J =1, demanded a specific contour for the Watson- 
Sommerfeld-Regge transformation, which could be translated by the need of 
subtraction in the dispersion relation [14]. That, in turn, could be represented 
by the derivative factor, {-k / 2)d / dins ^ in the odd case, Eq. (9) [31]. As we 
have shown, the DDR by Kang-Nicolescu and that by BKS have the same 
mathematical structure, once the cotangent operator (Section 3.1) could be 
well defined for the specific class of functions involved. We shall return to this 
important point in Section 5. 

At that time, the derivative approach was strongly criticized by several au- 
thors [17,18,19,20], in the sense that it had no practical interest, for instance, 
because "the mathematical condition for the convergence of the series excludes 
all cases of physical interest" [17]. However, it should be noted that some of 
these criticisms were connected with the physical situation at that time, some 
directed to attempts to extend the method to low energies (for example, as 
in [32]) and/or to some possible excessive optimism about the derivative ap- 
proach [33]. As effective contributions from all these criticisms we may mention 
the following results, which refer always to the even amplitude. 

Eichmann and Dronkers treated the case corresponding to a = 1, showing that 
Eq. (10) is valid only for some class of entire functions of Ins [17]. Specifically, 
with our notation of Sec. 3.1, they considered the series expansion of 
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and have proved that, if this series converges, then g{C,) must be an entire 
function of ^ and must satisfy \g{^)\ < Cexp{|^|}, C a constant. Conversely, 
if g{^) is an entire function of ^ and satisfies \g{^)\ < Cexp{(l — e)|^|}, e > 0, 
than the tangent series converges and may be represented by the integral 
relation (4), with sq = 0. Basically, the proof consisted in expanding g{C,) in 
Taylor series and compare the result with that obtained by performing the 
same expansion of g{^) in the integral (4). They have also shown that the 
method does not apply in the resonance region. 

Heidrich and Kazcs considered the parameter a and showed that the conver- 
gence of the series demands that a must lie between and 2 [18]. As before, 
it is assumed that ImF+(s)/s° may be expanded in Taylor series and that, 
in turn, integrated term by term. That leads to the evaluation of an integral 
depending on y = ^ — in the form 



-lncoth-|y| e(--^V|!~ + - / dy—--y^, 
TT 2 TT J smhv 

-oo ^ 

so that the first term is finite only for < a < 2. They have also shown that 
Eq. (6) is violated in the energy intervals between two branch points on the 
cut. It may be noted that a criterion of convergence based on the convergence 
radius of the Taylor expansion was later shown to be incorrect by Fischer and 
Kolaf [21]. 

Prom 1976 to 1987, Fischer and Kolaf, in a series of seminal papers [21], de- 
veloped a rigorous study of the derivative approach in connection with general 
principles and theorems from axiomatic quantum field theory [34]. The authors 
introduced classes of functions satisfying, among other properties, polynomial 
bound and the Froissart-Martin bound, and corresponding to a wider class 
of functions than that represented by entire functions in the logarithm of the 
energy. The analysis was focused on the region of asymptotic energies and the 
main formal result was that the derivative relations are valid if the tangent 
series is replaced by its first term. Although that demands the existence of 
the high-energy limits of certain physical quantities [21], the class of functions 
involved includes the majority of functions of interest in physical applications. 
For our purposes, we recall some results obtained by the authors about the 
role and range of the parameter a. For x = Ins, let F{x) be a function be- 
longing to the above mentioned class of functions [21]. For a = 1, the entire 
function J-'{z) which extends the function F{x) to the complex z plane, obeys 
the bound 
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|^(z)|<eel^l + C(e), 



for every e > 0, where C(e) is a constant that depends on e. On the other 
hand, for a 7^ 1, if the series 



tan 



2^" 



F{x) 



converges on some interval, then F[x) can be extended to an entire function 
for any a and in this case, the entire function obeys 



Therefore, the bound changes, unless x lies on the imaginary axis [21]. To our 
knowledge that is the only result in the literature that "quantify" the role 
of the parameter a. We shall return to this point in what follows. We also 
recall that differently from all the other works discussed here, Fischer and 
Kolaf make explicit reference to the fact that the tangent series represents the 
difference between lie F^/s and a constant associated with the subtraction. 
Moreover, they also discuss the contribution from the integral with zero and Sq 
as lower and upper limits, respectively (the high-energy approximation). They 
have also shown that for entire functions associated with fits, this contribution 
is divergent unless ImF(O) = (a result also presented in [18]). 

It should be also mentioned that Block and Cahn obtained a result for the 
DDR that does not depend on a, but under the assumption that the differential 
operator Z — d/dlns follows the condition \Z\ < 1 [2]. 

We conclude this discussion with the following comments. Concerning the 
range and validity of the parameter a, we have shown that the different state- 
ments, briefly quoted in Section 2, are not controversial since they are con- 
sequences of the different assumptions or classes of functions considered. As 
shown by Fischer and Kolar, for functions that arc entire in In s, the DDR are 
continuation of the IDR for all values of the parameter a. However, based on 
the above bounds, established for a = 1 and a ^ 1, we can not devise any 
practical advantage in using a 7^ 1, as stated before by Fischer and Kolaf. In 
this sense, we understand that the introduction of the parameter a by BKS 
is only an unnecessary complication. 

With the exception of the analysis by Fischer and Kolaf, neither of the works 
discussed in this section make any reference to the subtraction constant nor 
treat explicitly the contribution from the high-energy approximation repre- 
sented by the assumption sq ~* 0. In the following section we investigate the 
practical role of the subtraction constant and the influence of the lower limit 
of the IDR for a class of functions that are entire in Ins. 
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4 Tests with Pomeron-Reggeon pcirametrizations 

In order to quantitatively investigate the practical applicability of the DDR in 
substitution to IDR, we consider, as a framework, some standard parametriza- 
tions for the total cross section and present a detailed study on simultaneous 
fits to o'tot(s) and p{s). We start this Section with some comments about our 
choices concerning the parametrizations, ensemble of experimental data, en- 
ergy cutoffs and an outline on the strategies to be used in this study. We then 
present the fit results, followed by a discussion on the role of the subtrac- 
tion constant and of the high-energy approximation associated with the lower 
integral limit. 

4-1 Basic comments 

Presently, phenomenological analyses on the forward quantities atotis) and 
p{s), at high energies, are based on two kinds of dependences on the energy, 
namely power and power logarithmic functions and, therefore, entire functions 
in Ins. In the context of the Regge phenomenology they are associated with 
simple-pole Pomeron and secondary reggeons (s^''', < 7 < 1), double-pole 
(In s) and triple-pole (In^ s) contributions [16] . In order to treat in detail the 
applicability of the DDR and IDR by means of a complete example, we shall 
consider here, as a framework, only parametrizations based on simple-pole 
Pomeron and secondary reggeons, leaving the other cases for a forthcoming 
work. 

Although that Pomeron contribution, represented by the s"^ dependence, with 
e ~ 0.081 (Section 4.2), eventually violates the Froissart-Martin bound, we 
shall focus here only on its mathematical character, as representing a class 
of functions entire in Ins and, moreover, that can provide a good descrip- 
tion of the experimental data at the energies presently available. It should 
be also noted that other forms of DDR do not require the above bound, as 
demonstrated in [21]. 

Since all the Regge phenomenology is intended for the region of high energies 
[35] and, as mentioned previously, for particle-particle and particle-antiparticle 
(crossing) the pp and pp scattering correspond to the highest energy values 
with available data, we shall concentrate here only on these two processes. We 
shall return to this point in Section 4.4. We make use of the data sets on cxtot 
and p analyzed and compiled by the Particle Data Group [36]. The statistic 
and systematic errors have been added in quadrature. 

Since we are interested in the practical role of the subtraction constant as 
a free fit parameter, it is necessary to test both the different cutoffs in the 
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energy and the effect of the number of free fit parameters involved. To treat 
the former case, we shall consider in all the analysis two lower energy cuts: 
V^min ~ ^ GeV and v^min ~ GeV. In the later case, in order to get 
quantitative information on the sensitivity of the subtraction constant as a 
free fit parameter, we first consider an "economical" parametrization repre- 
sented by two reggeon exchanges (the Donnachie-Landshoff model) and after, 
an extended model with three reggeon exchanges. In the Regge context that 
means to consider degenerate and non- degenerate higher meson trajectories, 
respectively. 

Our aim in this section is to use these parametrizations in order to discuss the 
applicability of both the IDR, either with sq = 2m^ or sq = 0, and the DDR. 
Specifically, we want to establish how the results with the DDR deviate from 
those with the IDR and in which circumstances both approaches lead to the 
same results. Summarizing, our strategy shall be to consider simultaneous fits 
to (Ttot and p data from pp and pp scattering, with the following variants: (i) 
parametrizations with degenerate and non-degenerate trajectories; (ii) lower 
energy cuts at x/imin = 5 and x/Smin =10 GeV; (iii) subtraction constant: K 
— and X as a free fit parameter; (iv) DDR and IDR with both sq — 2vn? 
and with sq — Q. 

In what follows, we consider the approximation s = 4:{k^ + m'^) ~ 4A;^, so that 
the optical theorem, Eq. (2) reads 

(Ttot(s) = -■ (14) 



4-2 Degenerate meson trajectories 

The forward effective Regge amplitude introduced by Donnachie and Land- 
shoff (DL) has two contributions, one from a single Pomeron and the other 
from secondary Reggeons exchanges [37]. The parametrization assumes degen- 
eracies between the secondary reggeons, imposing a common intercept for the 
C — (02, 72) and the C — —1 {oj,p) and is given by 

aZ{s)=Xs^ + Ys~\ (15) 

aZ{s)^Xs^ + Zs-\ (16) 

Here, e = Q;p(0) — 1 and r] — 1 — aR(0), where Q;p(0) and q;r(0) are the 
Pomeron and Reggeon intercepts, respectively. From analysis of pp and pp in 
the interval 10 GeV < ^/s < 546 GeV, DL obtained the following values for 
the free parameters: X = 21.7 mb, Y = 56.08 mb, Z = 98.39 mb, e = 0.0808, 
rj = 0.4525. 
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In what follows, we first present the analytical connections between (Xtot and p 
with the DL parametrization, using both the IDR (with Sq = or Sq = 0) 
and the DDR and then, the results of global (simultaneous) fits to Utot and p, 
with K — or K as a, free fit parameter and energy cutoffs at 5 and 10 GeV. 



4-2.1 Analytical results 

The crossing even and odd amplitudes in Eq. (1) are obtained by using the 
parametrizations (15) and (16) and the optical theorem as given by Eq. (14). 
In this case, ImF+/_(s) oc s'^ , with — 1 < 7 < 1 and the integrals in (4) and 
(5) may be analytically evaluated. For sq = we obtain, 

P/ d.'^^^ = V-tan^. 

^ ^ 



For So = 2m^ fixed, the change of variables s' — y + sq and integration in y, 
lead to the result 



so 



where, from the above formulas, L{s) represents the contribution below the 
lower integration limit. 



^0 ,^ 



71 



and can be expressed either interms of hypergeometric functions or the corre- 
sponding series expansion in inverse powers of s [38] 



L{s) = ^ [2Fi(l, 7; 1 + 7; so/s) - 2Fi(l, 7; 1 + 7; -so/s)] 



With these expressions, and the DL parametrization, the connections between 
Ctot('S) and p{s) may be analytically determined. For sq = we obtain 



^±i^l^s-^cot L^)+XsHan(e'' 



p{s)aUs) ^ - 2 - -'-V2J V-2 
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-^^ -s 'tanir]— 

2 \ 2 



(17) 



and for fixed sq = 2m^, 



-^.-nan(,^) + iEf±^^^^ 
2 ^2^ TT^t^oV s(2j + 2-77) 

+ 27TTT^+ 2i + i- J UJ ' ^^^^ 

where the signs ± apply for pp (+) and pp (— ) scattering and, as before, the 
power series is associated with the contribution below the lower integration 
limit. 

The results with the DDR, as given by Eqs. (10) and (11) (or (12) and (13)), 

are exactly the same as that obtained by means of the IDR with Sq = 0, 
namely, Eq. (17). This is in agreement with our conclusion at the end of the 
Section 3.1: the only formal approximation involved in going from IDR to 
DDR is to take sq = 2m^ 0. 

Once we have established the analytical equivalence of the real parts obtained 
with both DDR and IDR for sq = 0, in what folllows we shall consider only 
the results provided by DDR, Eq. (17), and by IDR for sq = 2m^, Eq. (18). 



4.2.2 Global Fits 

We have developed global (simultaneous) fits to dtot and p with the DL param- 
etrization, using either the IDR with sq = 2m?, or the DDR, considering either 
X = or if as a free fit parameter and for energy cuts at 5 and 10 GeV. 
The fits have been performed with the program CERN Minuit and the errors 
in the fit parameters correspond to an increase of the by one unit. The 
numerical results and statistical information from all the four cases analyzed 
with energy cut at 5 GeV are displayed in Table 1, and the curves, together 
with the experimental data, are shown in Figs. 1 and 2, for X = and X as a 
free fit parameter, respectively. The corresponding results for the energy cut 
at 10 GeV arc presented in Table 2 and Figs. 3 and 4. We shall discuss these 
results in Section 4.4 together with those obtained by means of an extended 
parametrization, characterized by the non-degenerate trajectories. 
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Table 1 



Simultaneous fits to atot and p through the DL parametrization, \fs-^i^ = 5 GeV 
(238 data points), either with K as a free parameter or = and using IDR with 
lower limit sq = 2m^ and DDR. 





IDR with 


So = 2m? 


DDR 




K free 


K = 


K free 


K = 


X (mb) 


23.08 ± 0.28 


23.14 ± 0.28 


23.08 ± 0.28 


23.17 ± 0.28 


Y (mb) 


54.65 ± 1.3 


54.4 ± 1.3 


54.66 ± 1.3 


48.8 ± 1.0 


Z (mb) 


108.8 ± 3.5 


108.5 ± 3.5 


108.8 ± 3.6 


95.0 ± 2.4 


e 


0.0747 ± 0.0014 


0.0743 ± 0.0014 


0.0747 ± 0.0014 


0.0738 ±0.0014 


r] 


0.494 ± 0.010 


0.4946 ± 0.0099 


0.494 ± 0.010 


0.4673 ± 0.0079 


K 


17 ± 13 





179 ± 15 





XVF 


1.39 


1.39 


1.39 


2.01 


Table 2 










Simultaneous fits to crtot and p through the DL parametrization, i 


= 10 GeV 


(154 data 


points), cither with K as & free parameter or if = and 


using IDR with 


lower limit sq = 2m^ and DDR. 








IDR with 


So = 2m? 


DDR 




K free 


K = 


K free 


K = 


X (mb) 


21.62 ± 0.38 


21.78 ± 0.37 


21.62 ± 0.38 


21.83 ± 0.38 


Y (mb) 


65.7 ± 4.5 


59.5 ± 3.2 


65.7 ± 4.5 


50.6 ± 2.2 


Z (mb) 


113.4 ± 9.4 


102.5 ± 7.0 


113.4 ± 9.4 


85.7 ± 4.4 


e 


0.0816 ± 0.0018 


0.0806 ± 0.0018 


0.0816 ± 0.0018 


0.0800 ± 0.0018 




0.482 ± 0.019 


0.467 ± 0.017 


0.482 ± 0.019 


0.435 ± 0.014 


K 


116 ± 36 





287 ± 44 







1.17 


1.24 


1.17 


1.50 



4-3 Non-degenerate meson trajectories 



Analyses, treating global fits to a^ot and p, have indicated that the best re- 
sults are obtained with non-degenerate meson trajectories [39]. In this case 
the forward scattering amplitude is decomposed into three reggeon exchanges, 
F{s) — Fp(s) -I- Fa^if^{s) -\- tF^Ip{s), where the first term represents the ex- 
change of a single Pomeron, the other two the secondary Reggeons and r = +1 
(— 1) for ipp) amplitudes. Using the notation a]p(0) = 1 + e, «+(0) = 1 — ^7+ 
and q;_(0) = 1 — ?]_ for the intercepts of the Pomeron and the C = +1 and 
C — —1 trajectories, respectively, the total cross sections, Eq. (14), for pp and 
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T/s(GeV) T/s(GeV) 



Fig. 1. Simultaneous fit to atot and p through the DL parametrization, ^-^^^ = 5 
GeV, assuming if = and using either the IDR with sq = 2w? (sohd) or the DDR 
(dashed) . 




Vs(GeV) Vs(GeV) 



Fig. 2. Simultaneous fit to Ctot and p through the DL parametrization, \/Smiii ~ ^ 
GeV, with K as free parameter and using either the IDR with sq = 2to^ (solid) or 
the DDR (dashed): both curves coincide. 

pp interactions are written as 

atot(s) - Xs' + F+ 8-"+ + tY_ s-^- . (19) 

Since the formal structure of the parametrization is the same as in the case of 
degenerate trajectories, with the addition of one more power term, we shall not 
display the analytical formulas, but only the numerical and fit results, together 
with the experimental data. The numerical results and statistical information 
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100 




Vs(GeV) VsGeV 



Fig. 3. Simultaneous fit to cxtot and p through the DL parametrization, x/imjn = 10 
GeV, assuming K = {) and using either the IDR with sq = 1w? (soUd) or the DDR 
(dashed). 




from all the four cases analyzed with \fs^-^^ = 5 GeV are displayed in Table 
3 and the curves, together with the experimental data, are shown in Figs. 5 
and 6, for i^' = or as a free fit parameter, respectively. The corresponding 
results for \/injin = 10 presented in Table 4 and Figs. 7 and 8. 
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Table 3 

Simultaneous fits to atot and p through the extended parametrization, ^^^^ = 5 



GeV (238 data points), either with if as a free parameter ox K = 
with lower limit sq = 2m? and DDR. 


and using IDR 




IDR with 


So = 2m? 


DDR 




K free 


K = 


K free 


K = 


X (mb) 


19.61 ± 0.53 


20.12 ± 0.51 


19.61 ± 0.53 


21.44 ± 0.44 


y+ (mb) 


66.4 ± 2.0 


68.9 ± 2.1 


66.4 ± 2.0 


78.0 ± 2.0 


y_ (mb) 


-33.8 ± 1.9 


-32.4 ± 1.8 


-33.8 ± 1.9 


-38.8 ± 2.6 


e 


0.0895 ± 0.0025 


0.0874 ± 0.0024 


0.0895 ± 0.0025 


0.0814 ± 0.0021 


r}+ 


0.382 ± 0.015 


0.399 ± 0.014 


0.382 ± 0.015 


0.450 ± 0.012 


V- 


0.545 ± 0.013 


0.533 ± 0.012 


0.545 ± 0.013 


0.573 ± 0.015 


K 


-48 ± 15 





69 ± 18 







1.10 


1.14 


1.10 


1.64 


Table 4 

Simultaneous fits to atot and p through the extended parametrization, \/s^^^ = 10 
GeV (154 data points), either with K as a free parameter or K = and using IDR 
with lower limit sq = 2m^ and DDR. 




IDR with 


So = 2m^ 


DDR 




K free 


K = 


K free 


K = 


X (mb) 


19.57 ± 0.79 


19.70 ± 0.64 


19.58 ± 0.78 


21.088 ± 0.54 


Y+ (mb) 


66.0 ± 6.7 


67.4 ± 4.8 


66.0 ± 6.6 


86.850 ± 4.4 


y_ (mb) 


-29.2 ± 4.0 


-29.1 ± 3.9 


-29.2 ± 4.0 


-27.477 ± 6.0 


e 


0.0897 ± 0.0033 


0.0892 ± 0.0028 


0.0897 ± 0.0033 


0.0836 ± 0.0024 


V+ 


0.380 ± 0.033 


0.386 ± 0.024 


0.380 ± 0.033 


0.46256 ± 0.019 


V- 


0.520 ± 0.025 


0.519 ± 0.024 


0.520 ± 0.024 


0.50596 ± 0.040 


K 


-14 ± 48 





104 ± 58 





xVf 


1.10 


1.09 


1.10 


1.55 



4-4 Discussion 



Making use of power law parametrizations for the total cross section, we have 
presented the results of several simultaneous fits to o"tot(s) and p(s) from pp 
and pp scattering. Once established that the analytical results through IDR 
with So = and DDR are the same (Section 4.2.1), we have performed 16 
fits to the experimental data with the following variants: (1) degenerate and 
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Fig. 6. Simultaneous fit to (Jtot and p through the extended parametrization, 
V^min ~ ^ GcV, with K as free parameter and using either the IDR with sq = 2m^ 
(soUd) or the DDR (dashed): both curves coincide. 

non-degencratc trajectories; (2) v^j^in ~ ^ GeV and ^/Sjnin = 10 GeV; (3) 
IDR with So = and DDR; (4) K = and K as a free fit parameter. The 
results are displayed in Tables 1 to 4 and Figures 1 to 8. In what follows we 
focus the discussion on the effect of the high-energy approximation (sq — ^ 
in the IDR leading to the DDR) and on the role of the subtraction constant. 

For K = 0, the only difference between the IDR and the DDR approaches 
is associated with the high-energy approximation, analytically represented by 
integrations of Eqs. (4) and (5) from s' — to s' — sq — 2m^ f=:i 1.8 GeV. In 
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Vs(GeV) Vs(GeV) 



Fig. 7. Simultaneous fit to atot and p through the extended parametrization, 
V^min = 10 GeV, assuming K = Q and using either the IDR with sq = 2m^ (soUd) 
or the DDR (dashed). 




Vs(GeV) Vs(GeV) 



Fig. 8. Simultaneous fit to crtot and p through the extended parametrization, 
V^min ~ 10 GcV, with K as free parameter and using either the IDR with sq = 2m? 
(solid) or the DDR (dashed): both curves coincide. 

principle, it might be expected that this contribution could affect the fit results 
only in the low-energy region, roughly below ^ 50 GeV. However, the fit 
procedure is characterized by strong correlations among the free parameters 
and, therefore, the contributions from any region may be "communicated" to 
other regions and that is exactly what our results show. In fact, for K — 
0, independently of the parametrization and the lower energy cut used, the 
leading contribution at high energies (represented by the value of the Pomeron 
intercept, e) is different when using IDR or DDR. The above approximation 
(DDR) results in a value for the intercept that is smaller than that obtained 
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with the IDR. Specifically, in the case of the DL parametrization (Tables 1 and 
2). the reduction reads 0.67 % for v^j^i^ = 5 GeV and 0.7 % for ^/Sj^jn = 10 
GeV (see crtot(s) in Figs 1 and 3). The effect is more significant in the case of 
the extended parametrization (Tables 3 and 4), with a reduction of 6.9 % for 
v^^i^ = 5 GeV and 6.3 % for y^^.^ = 10 GeV (see (Jtot(s) in Figs. 5 and 7). 
That is an important and novel result, which shows that, in practice, for the 
simple-pole Pomeron and secondary rcggcons, the high-energy approximation 
enclosed in the DDR affects the fits results even at the asymptotic energies. 
As mentioned, that is a consequence of the fit procedure and the correlations 
among the free parameters. 

Now, the practical role of the subtraction constant may be investigated by 
comparing the results obtained with K = and those with as a free fit 
parameter. The striking result here is the fact that, independently of the 
parametrization and the energy cut used, the numerical results obtained for 
the fit parameters with the IDR and with the DDR are the same up to 3 
significant figures, including the 'X^ /F. For example, from Tables 1 and 2, for 
the case of the DL parametrization, we have obtained in both cases (IDR with 
So = and DDR) the same values: e = 0.0747 ± 0.0014, / F = 1.39 for 
V^^i„ = 5 GeV and e = 0.0816 ± 0.0018, / F = 1.17 for V^min = 10 GeV 
(sec Figs 2 and 4, respectively). For the cases of the extended parametriza- 
tion. Tables 3 and 4, we have obtained e = 0.0895 ± 0.0025, x^ I F = 1.10 for 
Vs^in = 5 GeV and e = 0.0897 ± 0.0033, x^/F = 1.10 for y^^;^ = 10 GeV 
(see Figs 6 and 8, respectively). 

All these numerical results can be understood in an analytical context if we 
investigate the correlation between the subtraction constant and the contribu- 
tion below the lower integration hmit (referred to in Sec. 4.2.1). To this end, 
let us return to the product of p{s) by atdis) in the case of IDR with sq — 2m^, 
as given by Eq. (18), or the corresponding formula with non-degenerate trajec- 
tories. In both cases the terms with entire inverse powers of s can be expressed 

by 

[K + A]- + 0(l/s% 

s 

where K is the subtraction constant and A comes from the series expansion, 
which is associated with the integration from s' = to s' = 2m^ (Sec. 4.2.1). 
Therefore, at large enough energies, this contribution below the lower integra- 
tion limit is expected to be absorbed in the subtraction constant. In fact, with 
degenerate (de) and non-degenerate (non-de) trajectories, these contributions 
read 

1 f2Xso^ {Y + Z)so'-^ \ 

^de = - — r- \ z , 

n \ 1 + e 1-7] J 



24 



A. 



non— de 



1 / 2Xso^+' ^ 2Y+so'-^+ 
7r\l + e l-?7+ 



) 



and may be estimated with the central values of the fit parameters from Tables 



Ade ~ 162 = 5 GeV), Aje ~ 170 {^^^ = 10 GeV), 

Anon-de ~ US (v^^;, = 5 GeV), Ano„_de ~ 117 (v/i^i„ = 10 GeV), 



From those Tables, we see that these values are in complete agreement with the 
differences between the fitted values of the subtraction constants determined 
with DDR and IDR for sq = 2m^. Therefore, we conclude that in the case 
of DDR we have an "effective" subtraction constant which can compensate 
(analytically and numerically) the effect of the high-energy approximation. 

That is another important and novel result, showing explictly the practical 
role of this parameter: for simplc-polc pomeron and secondary reggeons once 
the subtraction constant is considered as a free fit parameter, the DDR are 
completely equivalent to the IDR with fixed (non-zero) lower limit. We also 
note that our qualitative and quantitative conclusions about the effects of the 
high-energy approximation and the subtraction constant are the same with 
both the degenerate and non-degenerate trajectories and both energy cuts, 5 
and 10 GeV. 

To conclude this Section, let us call attention to some critical aspects related 
with the use of the DL and the extended parametrizations. From a statistical 
point of view, the extended paramctrization leads to a better description of 
the experimental data (Tables 1-4). Also, the p data at 546 GeV and 1.8 
TeV are well described with the extended paramctrization, but not in the DL 
case. The essential difference between these two parametrizations concerns the 
additional secondary reggeon, represented by the dependence s"''', < 7 < 
1, which goes to zero as the energy increases. We conclude that the splitting 
of the trajectories allows the free parameters from the secondary reggeons 
to fit the data at low energies (between 5 GeV and, let us say, 100 GeV), 
giving more freedom for the parameters associated with the Pomeron to fit 
the data at the highest energies. That seems to be well known and also well 
accepted in the literature [39]. However, we can not forget that the Rcgge 
phenomenology is intended only for asymptotic energies. Specifically, the basic 
contribution comes from the asymptotic form of the Legendre Polynomial, 
namely Pi{x) —>■ x'' as x — > 00, where x is translated to the energy s through 
the crossing symmetry and / to the trajectory of the Pomeron or Reggeon [35]. 
On the other hand, from the above discussion, the secondary reggeons act in 
a region where the total cross section reaches its minimum (followed by an 
increase) and the p parameter cross the zero (becoming positive). Even if we, 
optimistically, concentrate in the above region, ^/s = 5 to 100 GeV, that 
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region, certainly, has nothing to do with an asymptotic concept. Based on 
these facts, wc understand that, although the secondary reggeons can be used 
as suitable parametrizations on statistical or strictly fit grounds, one must 
be careful in attempting to extract well founded physical results from these 
asymptotic forms used at so low values of the energy. On the other hand, that 
is not the case for the subtraction constant as a free fit parameter, since it has 
well defined mathematical bases, which are related with polynomial bounds. 



5 Conclusions and final remarks 

In this work we have presented a review on the different results and statements 
in the literature related with the replacement of IDR by DDR, and a discussion 
connecting these different aspects with the corresponding assumptions and 
classes of functions considered in each case. 

By means of a formal and analytical approach, we have demonstrated that 
the subtraction constant is preserved when the IDR arc replaced by the DDR 
and that, for the class of functions entire in his, the DDR do not depend on 
any additional free parameter (except for the subtraction constant). We have 
stressed that the only approximation involved in this replacement concerns the 
lower limit sq in the IDR: the high-energy condition is reached by assuming 
that So = 2m^ 0. 

We have investigated the practical apphcability of the DDR and IDR in the 
context of the Pomeron-reggeon parametrizations, with both degenerate and 
non-degenerate higher meson trajectories. By means of global fits to atotis) 
and p(s) data from pp and pp scattering, we have tested all the important 
variants that could affect the fit results, namely the number of secondary 
reggeons, energy cutoff, effects of the high-energy approximation and the sub- 
traction constant and the derivative approach with DDR and IDR with fixed 
Sq. Our results lead to the conclusion that the high-energy approximation and 
the subtraction constant affect the fit results at both low and high energies. 
This effect is a consequence of the fit procedure, associated with the strong 
correlation among the free parameters. 

A striking novel result concerns the practical role of the subtraction constant. 
We have shown that, with the Pomeron-reggeon parametrization, once the 
subtraction constant is used as a free fit parameter, the results obtained with 
the DDR and with the IDR (with finite lower limit, Sq = 2m^) are the same up 
to 3 significant figures in the fit parameters and / F- Analytically this effect 
is due to the absorption of the high-energy approximation in an "effective" 
subtraction constant. This conclusion, as we have shown, is independent of 
the number of secondary reggeons (DL or extended parametrization) or the 
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energy cutoff {^/s = 5 or 10 GeV). 



We fiave called the attention to the fact that the subtraction constant has 
well founded mathematical bases, since it is a consequence of the polynomial 
bounds in the scattering amplitude. On the other hand, the use of the asymp- 
totic Regge forms for detailed fits at finite energies (minimum on o'tot(s)) is not 
formally justified. We understand that our results suggest that, presently, it is 
very important to look for usefuU and well founded theoretical results at finite 
energies and that one must be careful in attempting to predict asymptotic 
behaviors based on asymptotic formalisms applied at finite energies. 

In order to treat a complete example, with the 16 variants referred to in 
the text, we have considered here only a simple-pole representation for the 
Pomeron. Although belonging to the class of functions entire in In s, this 
choice is outside the more general class of functions considered by Fischer 
and P. Kolaf, because it violates the Froissart-Martin bound at the asymp- 
totic energies. We are presently investigating other possibilities, as dipole and 
tripole contributions. Anyway, there is an interesting aspect concerning these 
forms, which we shall mention here. The possibility that the odd amplitude 
had a Ins contribution (Odderon) has been investigated in Refs. [14,26], by 
means of the DDR. However, this contribution in the IDR, Eq. (5), with Sq = 
is divergent. On the other hand, that is not the case if we use the DDR as 
expressed by Eqs. (9) or (13). Since we have shown that the only approxima- 
tion in going from IDR to DDR is to take sq ~^ in the integral relation, the 
above results seems inconsistent. We think that this may be associated with 
the formal substitution of the tangent by the cotangent operator, as referred 
to in Sec. 3.1. We are presently investigating this subject. 

The DDR introduced by Bronzan, Kane, and Sukhatme [13], and later gener- 
alized for an arbitrary number of subtractions by Menon, Motter and Pimentel 
[23] bring enclosed the parameter a, which we have shown does not need to 
take part in the calculation. We understand that this dependence is not nec- 
essarily wrong if, in practical uses, one takes a — 1. Therefore, the presence of 
this parameter is only an unnecessary complication. Certainly, considering a 
as a free parameter at finite energies may improve the fit results [13,24], but 
it seems to us difficult to justify its mathematical and, mainly, its physical 
meaning. 

Note added in proof. After the submission of this work, the effect of the absorp- 
tion of the high-energy approximation (sq — 0) by the subtraction constant 
has been confirmed by other authors [40]. 
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